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In the past decade, several neutron reflectometry methods for determining the modules and phase 
of complex reflection coefficient of an unknown multilayer thin film have been worked out among 
which the method of variation of surroundings and reference layers are of highest interest. These 
methods were later modified for measurement of polarization of reflected beam instead of the mea- 
surement of intensities. In their new architecture, these methods not only suffered from the necessity 
of change of experimental setup, but also another difficulty was added to their experimental im- 
plementations which was related to the limitations of the technology of the neutron reflectometers 
that could only measure the polarization of the reflected neutrons in the same direction as the 
polarization of the incident beam. As the instruments are limited, the theory has to be optimized 
so that the experiment could be performed. In a recent work, we have developed the method of 
variation of surroundings for one directional polarization analysis. In this new work, the method of 
reference layer with polarization analysis has been optimized to determine the phase and modules of 
the unknown film with measurement of the polarization of reflected neutrons in the same direction 
as the polarization of the incident beam. 

PACS numbers: 61. 12. Ha Neutron reflectometry 



I. INTRODUCTION 

As an atomic scale probe, neutron reflectometry with 
polarized neutrons has a vast application to the study of 
layered media and thin films. Each material has a unique 
scattering length density (SLD) for neutrons which can 
be used as a unique characteristic for determining the 
type and thickness of that material 1 . In other words, re- 
trieving the nuclear scattering potential of neutrons from 
reflectivity data can be used as a powerful tool for iden- 
tifying the unknown samples. However, measuring the 
intensity of the reflected neutrons from the sample would 
not lead to a unique result for the SLD of the sample. In 
the scattering process, the reflectivity R(q) =| r{q) | 2 is 
measured and the information of the phase of the com- 
plex reflection coefficient r(q) is lost. As a result, two 
different layer with different phase information can have 
the same reflectivity. Hence, the full knowledge of the 
phase of complex reflection coefficient is necessary in or- 
der to retrieve a unique result from the reflectometry 
experiments 1,2 . 

In the paste decades, several theoretical methods for 
resolving this so-called phase problem have been worked 
out such as variation of surroundings medium 3,4 which 
makes use of the controlled variations of scattering length 
density of the incident and/or substrate medium or the 
method of reference layers 5 ' 6 which is based on the in- 
terference between the reflections of a known reference 
layer and the unknown surface profile. Either of the 
methods suffer from the experimental difficulties regard- 
ing the change of the substrate or reference layers for 
each reflectivity measurement 4 ' ' . This inefficiency was 
recovered by using polarized neutrons and magnetic lay- 
ers as substrate 8,9 or reference layer 10-12 . Later on, these 
methods were enhanced to measure the polarization of 
the reflected beam to determine the phase of reflection 13 . 



Although the polarization based approaches had been 
theoretically proven to be efficient, their experimental 
implementation was limited by the ability of reflectome- 
ters in the polarization measurement direction. As these 
methods required at least two measurements of polar- 
ization of the reflected bean in different directions 1 1-13 , 
their experimental implementations were not practical 
with reflectometers which can only measure the polar- 
ization of the reflected neutrons in the same direction as 
the incident beam. 

In a recent work 14 , we developed the method of vari- 
ation of surroundings with one directional polarization 
analysis and discussed about possibility of experimental 
implementation of the method. In this new worked, we 
have optimized the method of reference layers based on 
the measurement of the polarization of reflected beam in 
the same direction as the incident neutrons. 

In section II, the foundation of the optimized reference 
layer is formulated. Section III, deals with numerical ex- 
amination of the method for an unknown sample. The 
SLD of the sample is also retrieved from the phase in- 
formation which was obtained from optimized reference 
layer method. The paper ends with further discussions 
about the experimental challenges of the method. 



II. OPTIMIZED REFERENCE LAYER METHOD 

Scattering of neutrons from a sample is described by 
optical potential, v{z) — 2irh 2 p(x)/m, where p(x) = 
p n ± p m is the scattering length density of the sample 
as a function of its depth. p n is the nuclear part of the 
SLD and p m = (m/2irh 2 )pB is the magnetic part. The 
Plus (minus) signs denotes the incident beam polarized 
parallel (anti parallel) to the local magnetization. The 
scattering of neutrons from the sample is described by 
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FIG. 1. Pictorial demonstration of the arrangement of the 
sample. An unknown bilayer film is mounted on top of the 
known magnetic reference layer which is deposited over an 
infinite substrate 



the one dimensional Schrdinger equation: 

{d 2 x + (q 2 -47rp(x))}iP( qi x) = 



(1) 



where q is the incident neutron wave number in the x 
direction. The reflection (r±) and transmission (t±) co- 
efficients of the sample can be determined from Eq. 1 
using the transfer matrix method 2 : 
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where L is the thickness of the sample and j = (1 — 
Ampj/q 2 ) 1 / 2 with j = f,h is the refractive index of 
the fronting and backing medium respectively (pf and 
Ph are the SLDs of the fronting and backing medium). 
(A, . . . , D) are the elements of the transfer matrix which 
are uniquely determined as a function of the SLD of the 
sample. The reflection coefficient of the sample in terms 
of the transfer matrix elements is also written as follows: 
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The reflectivity depends on a± , f3± and 7± in terms of a 
new quantity £±: 
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Suppose the sample is composed of two part, an un- 
known part which is mounted on top of the known mag- 
netic reference layer which is magnetized in the -In- 
direction (Fig. I). The transfer matrix of the whole sam- 
ple is written as the multiplication of the transfer matrix 



where (a ... d) is for the unknown part of the sample and 
(w ... z) belong to the known part. Using Eq. 5,6 , X± 
is modified as follows: 
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where the subscript k (u) denotes the known (unknown) 
part of the sample and the tilde represents the mirror- 
reversed sample which is obtained by the interchange of 
the diagonal elements of the transfer matrix, (A «-» D). 
The superscript // denotes a sample with the same 
medium at both sides with refractive index /. 

The polarization of the reflected beam, (j? Xl p y ,p z ), is 
written in terms of the polarization of the incident beam, 
(p£,?M) 10 ' 15 : 



2(rf )*r f _ h (p° x + ip°) 
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Up to now, we have defined all the preliminary parame- 
ters. Here, we are going to perform two different polar- 
ization measurements which satisfy our goals: 

First measurement: We suppose the incident neu- 
trons to be fully polarized in the y direction. The polar- 
ization of the reflected beam, Pi , is also measured in the 
same direction as the incident beam. From Eq. 8 for the 
measured polarization, we have: 



Pi = 1- 
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which is completely determined from the knowledge of 
the known part of the sample. 

Second measurement: For the second case, we con- 
sider a non-polarized incident beam, rotate the sample 
by 90° in the y — z plane and measure the polarization of 
the reflected neutrons in the the y direction. This mea- 
surement is the same as measuring the polarization in 
the direction of the local magnetization of the reference 
layer, z, for the non-rotated sample. This choice for the 
polarization analysis direction would not contradict with 
our first assumption of fixed experimental setup. From 
Eq. 8, the second polarization can be written as: 
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As the knowledge of Pi and Pi are experimentally 
known, is determined from the following quadratic 
equation: 
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This quadratic equation has two different solutions; the 
physical solution is selected using the fact that E + > 2. 
By knowing e£ + v '" 



the E_ is also determined as follows: 



I -Pi 



(13) 



Using Eq. 7 and the knowledge of Ej' 1 



j ± , the parameters 
of the unknown part of the sample are related to each 
other as follows: 



TABLE I. Information of the sample 



Layer Thickness (nm) p n (10 4 nm 2 ) p m (10 4 nm 



NiO 20 8.84 

Fe 2 3 15 7.26 

Co 20 2.23 

Si Semi Infinite 2.08 
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By using Eq. 14 and the fact that 7 = \JafB — l 6 , the 
jff parameter is determined from the following quadratic 
equation: 

(aia 2 - 1)7 2 + (aib 2 + a 2 &i)7 + {b\b 2 - 1) = (16) 

The complex reflection coefficient is then determined as 
a function of : 
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The quadratic equation 16 has two different solutions 
which only one of them is physical. The physical so- 
lution can not directly be determined from Eq. 16. In 
order to determine the physical solution, we put the two 
solutions into Eq. 17 and select the physical one from the 
behavior of the reflection coefficient: 



< Re{r) < 1, < Im(r) < 1 
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As both of the two solutions satisfy Eq. 18, the physi- 
cal solution should be determined from their behavior in 
q — » (Eq. 19). When the physical solution in q — > is 
determined, the solution in larger q values is also deter- 
mined based on the continuity of the reflection coefficient 
along the whole range of q values. 
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FIG. 2. Polarizations of the reflected neutrons for the two 
different measurements 



III. NUMERICAL EXAMPLE 

In this section, we numerically examine the optimized 
reference layer method to certify the reliability of the 
method. We consider (NiO / Fe 2 Os) bilayer as the un- 
known sample which is mounted on top of a Cobalt layer 
as magnetic reference. The whole sample is deposited 
over a silicon substrate. The SLD and thickness of the 
layers are listed in table. I. According to the the opti- 
mized reference layer method, for the first step, a y- 
polarized neutron beam, P — (0,1,0), was radiated to 
the sample and the polarization of the reflected beam, Pi , 
was measured in the same direction. For the second mea- 
surement, first the sample is rotated by 90° in the y — z 
plane and then by using a non-polarized incident beam, 
P = (0,0,0), the polarization of the reflected beam, P 2 , 
is measured in the z direction. Using the measured po- 
larizations in Eqs. 9, 11 and following the Eqs. 12-17, the 
complex reflection coefficient of the mirror image of the 
unknown part of the sample while it is surrounded by 
vacuum at both sides, is determined. 

Fig. 2 shows the polarization of the reflected beam for 
the two different measurements. Fig. 3 also demonstrate 
the imaginary part of the complex reflection coefficient 
of the sample which is retrieved from Eq. 17. As it is 
illustrated in the figure, the physical solution alternates 
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FIG. 3. Imaginary part of the complex reflection coefficient 
which is determined from Eq.f7. The physical solution al- 
ternates between the two different solutions of the quadratic 
equation (blue and red curves). The solid black line shows 
the physical solution 




FIG. 4. Reflectivity of the mirror reversed sample with vac- 
uum surroundings. The retrieved reflectivity (red curve), 
clearly corresponds to the measured reflectivity 



between these two sets. The physical solution is selected 
based on the continuity of the data and the fact that as 
q — > 0, the imaginary part of the complex reflection co- 
efficients, Im(r) — > from the negative 4 . The retrieved 
reflectivity of the sample is also illustrated in Fig. 4 which 
truly corresponds to the measured reflectivity from the 
experiment. 

Knowing the information of the complex phase of re- 
flection, the SLD of the sample can be determined by 
solving the inverse Schrdinger equation for the scatter- 
ing potential. In neutron reffectometry problems, the 



FIG. 5. Retrieved Optical Potential v of the sample from the 
data of Eq.17 and using the Sacks Code 1 ' . 



SLD of the sample is retrieved by solving the GclfandLc- 
vitanMarchcnko integral equations 2 ' 3 ' 16 . Sacks et al. 1 ' 
also developed a numerical code in which the real and 
imaginary part of the complex reflection coefficient are 
used as input and the scattering optical potential is de- 
livered as output. Fig. 5 shows the retrieved SLD of our 
sample which is determined by using the Sacks code. One 
can see that the retrieved SLD corresponds to the mirror 
image of the unknown part of the sample with vacuum 
surroundings at both sides. 



IV. CONCLUSION 

Available neutron refiectometers are limited in the 
measurement direction of the polarization of the reflected 
beam according to which the polarization of the reflected 
beam can only be measured in the same direction as 
the polarization of the incident beam. The polarization 
based theoretical approaches which deal with the deter- 
mination of the missing phase of reflection, suffer from 
this experimental deficiency. In this paper, we have op- 
timized the method of reference layers based on the one 
directional polarization analysis and resolved this prob- 
lem. 

The optimized reference layer method deals with two 
different polarization measurements in the same direction 
as the polarization of the incident beam. Moreover, the 
method resolves the the necessity of the change of exper- 
imental setup for polarization measurements which is an 
important issue in experiments. In the optimized refer- 
ence layer method, the information of the complex reflec- 
tion coefficient is retrieved only with polarization anal- 
ysis. And in contrast to most of other method, there is 
no need to the measurement of reflectivity. The method 
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provides us with the modules and phase of the complex 
reflection coefficient and results to the unique determina- 
tion of the scattering length density of the mirror image 
of the unknown sample. 

In order to check the reliability of the method, a nu- 



merical simulation was performed and the knowledge of 
the missing phase of the reflection and as a consequence, 
the SLD of the sample was truly retrieved. The numerical 
results certify that the method is reliable at this level and 
we look forward to its experimental implementations. 
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